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Abstract 

The quasiparticle time blocking approximation (QTBA) is applied to describe El 
excitations in the even-even tin isotopes. Within the model pairing correlations, two- 
quasiparticle (2q), and 2q(8)phonon configurations are included. Thus the QTBA is an 
extension of the quasiparticle random phase approximation to include quasiparticle- 
phonon coupling. Calculational formulas are presented in case of neutral excitations in 
the spherically symmetric system. The main equations are written in the coordinate 
representation that allows to take into account single-particle continuum completely. 
The El photoabsorption cross sections have been calculated in nuclei n6,i20,i24gj^_ 
has been obtained that the 2q(8)phonon configurations provide noticeable fragmentation 
of the giant dipole resonance resulting in appearance of significant spreading width. 
The results are compared with available experimental data. 
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1 INTRODUCTION 


Theoretical description of giant mnltipole resonances (GMRs) and resonance structures in 
magic and open-shell nuclei has a long history and remains a problem of current importance. 
The main tools in solving this problem within the framework of microscopic approach are the 
random phase approximation (RPA) and the quasiparticle RPA (QRPA) which is straight¬ 
forward generalization of the RPA to include pairing correlations. However, despite these 
models enable one to describe mean energies and total strengths of the GMRs, they fail to 
reproduce the total widths of the resonances and their hne structure (see, e.g.. Ref. PP). 
One of the reasons is that RPA and QRPA do not provide a mechanism producing spreading 
width pJ^ which gives a considerable contribution in the total widths of the GMRs. 

Simplest mechanism of this type is the coupling of the quasiparticles to phonons being 
superpositions of the one-particle-one-hole (Iplh) or the two-quasiparticle (2q) configura¬ 
tions. As applied to structure of the even-even nuclei, the concept of the quasiparticle-phonon 
coupling (QPG, see Ref. |2j) enables one to take into account lplh®phonon and 2q®phonon 
conhgurations in addition to Iplh and 2q ones incorporated within the RPA and the QRPA. 
A series of models has been developed to go beyond the RPA and the QRPA by means of 
inclusion of this mechanism (see Refs. Q, Q-ISI and references therein). Recently, new 
model has been formulated (see Ref. m in which pairing correlations, 2q, and 2q®phonon 
configurations are included. This model is based on the Green function (GF) formalism 
that mainly determines its name: quasiparticle time blocking approximation (QTBA). The 
QTBA is a generalization of the method of chronological decoupling of diagrams (MGDD) 
developed in Ref. HOI to describe excited states of the even-even nuclei without pairing. 

The hrst aim of this paper is to present calculational formulas obtained from the gen¬ 
eral ones of the model [Oj making use of certain approximations in case of neutral excitations 
in the spherically symmetric system. In particular, in the formulas presented zero-range 
forces are adopted as an effective interaction and the Bardeen-Gooper-Schrieffer (BGS) ap¬ 
proximation is used to determine quasiparticle energies and wave functions. The equations 
obtained are a system of coupled equations corresponding to the excitations in particle-hole, 
particle-particle, and hole-hole channels. The basic equations are written in the coordinate 
representation that allows to take into account single-particle continuum completely. Notice 
that inclusion of the continuum is of particular importance to describe correctly the widths 
of the GMRs. 
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Our second aim is to test the QTBA in calculations of the electric dipole excitations 
in the nuclei with pairing. With that end in view we have chosen tin isotopes ii6,i2o,i24gj^ 
We present El photoabsorption cross sections calculated within QRPA and QTBA. Calcu¬ 
lated integral characteristics of the giant dipole resonance (GDR) are compared with the 
experimental data. 

The paper is organized as follows. In Sec. |21 general scheme of the approach is pre¬ 
sented. In Sec. Olthe basic equations of the QTBA are transformed to channel form in the 
coordinate representation. In Sec. IH these equations are further reduced to the equations for 
partial components. Formulas for the correlated propagator of the QTBA in the coordinate 
representation are presented in Sec. El In Sec. El the equations for the matrix elements of 
the amplitude of quasiparticle-phonon interaction are obtained. These quantities enter the 
equations for the correlated propagator of the QTBA in terms of the reduced matrix elements 
which are presented in Sec. 13 In Sec. IHl we describe our calculational scheme and present 
the results obtained for the photoabsorption cross sections in three tin isotopes. Conclusions 
are drawn in the last section. 

2 GENERAL FRAMEWORK 

The basic quantity, which determines the physical observables in the QTBA, is the nuclear 
polarizability n(a;). More precisely, it determines distribution of the transition strength 
caused by an external held V^. In the representation of the single-quasiparticle basis func¬ 
tions, which will be specihed below, n(a;) is dehned as 

nn = -5 E . (1) 

^ 1234 

where e is the effective charge operator, a; is complex energy variable, is the effective 

response function. is a solution of the following Bethe-Salpeter equation 

-^12,34(^) ~ ^12,34 (^) ~ ^12,56 (^) -^56,78 -^78,34(^) ’ (2) 

5678 

where A (in) is a correlated propagator, JF is an amplitude of the effective interaction. In 
particular, the strength function S(B) which is frequently used for the description of nuclear 
excitations is expressed in terms of the polarizability as 

A(E) =--Imn(E + iA), (3) 

77 
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where A is a smearing parameter. 

The Eq. (j21) for the effective response function is quite general. It was obtained 
in Ref. jHj for the Fermi systems with pairing within the framework of generalized Green 
function formalism in which the normal and the anomalous GFs are treated in a unihed way 
in terms of the components of generalized GFs in a doubled space. The physical content 
of the model is determined by the choice of the propagator A{uj). In particular, the QRPA 
equation for the has the same form m with uncorrelated propagator taken instead 

of A{uj). The equations dehning A{uj) within the QTBA are given in Ref. [H]. In the 
present paper we consider the version of QTBA in which the correlated propagator includes 
contributions of the 2q and 2q(8)phonon conhgurations. 

As a hrst step we have to determine the single-quasiparticle basis functions 
which form matrix representation of the theory. For the Fermi systems with pairing these 
functions are dehned in the doubled space spanned by the coordinates y = {x,x}, where 
symbol x = {r, a, r} includes the spatial coordinate r, the spin a, and the isospin r variables, 
X = ±1 is an additional index introduced for denoting the different components of the single¬ 
quasiparticle functions in the doubled coordinate space. Index 1 = {A 1,771} of the doubled 
conhguration space includes index Ai of the usual single-particle configuration space and the 
index 771 = ±1 which is the sign of the eigenvalue corresponding to "01 (7/). Namely, we have: 
7 -f-01 = Fi -01 where Ti is single-quasiparticle Hamiltonian, Ei = 771 E\^ = |Fi|. In case 

of the spherically symmetric system, we are interested in, the index A can be represented 
by the following set: A = {(A),r?7.} where (A) = {ta, tt,,/, j}, m is the projection of the total 
angular momentum, and we have: Ex = E(^xy 

We will use the approximation corresponding to so-called canonical basis representa¬ 
tion of the functions 'ipihj) (see Ref. jTT] for details). To determine functions 'ipiiy) within 
this approximation let us note that in the matrix form the Hamiltonian Ti reads 


n={^ ^ ^ \ 

I -A* y-h* j 


( 4 ) 


where h = h{x,x') is the single-particle Hamiltonian, A = A{x,x') is the operator of the 
pairing held, p is the chemical potential. Let {v 3 a(^)} complete set of orthonormal 

eigenfunctions of the Hamiltonian h{x,x'): h(p^ = We will assume that the op¬ 
erator A{x,x') has the canonical form in the same basis corresponds to the 

state-dependent version of the BGS approximation (see, e.g.. Ref. 113). or to the so-called 
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approximation of the diagonal pairing. In this case for the spherically symmetric system we 
have: 

^A,+ (a;, +) = Ma ) V’A,+ (a;, -) = (-1)^+^'+”^ n;, (pi(x) , 

^/’A,-(a;,-) =U),ipl{x), ?/;A,-((r,+) = , 

where A = {(A),—m}, v\ and u\ are real numbers which satisfy the following conditions: 
“a = VT^ ^ 0, vl = vl^ I [see Eqs. for an explicit definition of v\ and 

u\ within the BCS approximation]. The choice of the phase factors is determined by the 
formulas: 



= (-1)'+J+”+U-4,!(j;), (8) 

where n = r/r, x = {r, —a, r}. 

In applications of the theory it is convenient to use another basis functions which 
differ from the functions by a unitary transformation. Let us introduce a matrix 


0l2 = 0 




= h 


ni,V2 ,+i 




(9) 


This matrix is real and orthogonal, and consequently it is unitary. So wave functions 'ipi{y) 
defined through the single-quasiparticle basis functions 'ipi{y) by the formula 


i’liy) = '^Oi2'il^2{y) ( 10 ) 

2 

also form a complete set of the orthonormal functions. We will use just the set {'ipiiy)} as 
the set of basis functions. This does not lead to an inconsistency since the single-particle GF 
G(£) = (e — 'H)~^ is diagonal both in {V’i(2/)} and in {'0i(2/)} representation, and hence the 
formulas for the correlated propagator (see Ref. jH]) are the same in both representations. 

To describe dynamics of the system and to calculate the polarizability and the strength 
function we start from the equation for linear response matrix (LRM) A. Notice that in the 
present work the term LRM is used instead of the frequently used one “density matrix 
variation in an external field” because it is more correct in our notations. In the coordinate 
representation the equation for A is obtained by the convolution of the equation for effective 
response function with an operator of the renormalized external field eV^: 

R{yi,y2]^) = A^iyi,y2; uj) - j dy3dy4dy5dyeA{yi,y2; y3,y4; u) 

X -^(i/ 3,2/4; 2/5,1/6) A(|/5,2/6; a;), (11) 
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where 


( 12 ) 


Hvi, ‘ 1 ') = - i’l (yi) i’2{v2) (cv °)„, 

1234 


A°(i/i, 2/2; o;) = - j dx:idx^[A{yi,y2] X^+,X2^+] Uj) - A{yi,y2] Uj)] 

X V °{ X 4 , X 3 ). (13) 

It is assumed that the correlated propagator of the model A{(jj) is initially calculated in 
conhguration space and then is transformed to coordinate representation: 


A(2/i,2/2; 2/3,2/4; = X]^i*(2/i) 4(2/2) ^3(2/3) ^4(2/4) ^12,34(0;) . 

1234 

Components of the external field are 

V°(xi,X2) = V^{xi+,X2+) = -V^{x2-,xi-) = ^4(a^i+)4(3^2+) (eC°)i2 , 


(14) 


(15) 
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and it is supposed that C°(xi+,a; 2 —) = C°(xi—,X 2 +) = 0. 

In terms of the LRM the Eq. 0 for the polarizability reads: 

11 ( 0 ;) = / dxi dx 2 {x 2 , xi) A(xi+, X 2 +; a;). 


( 16 ) 


To determine general form of the effective interaction in our approach let us note 
that within a self-consistent scheme the amplitude JF in Eq. (HH) can be defined as a second 
order functional derivative of some energy density functional S[TZ] where TZ is the extended 
density matrix (7^(2/l,2/2) = E 4,-i see, e.g., Ref. [H]): 


-^(221,2/2; 223,2/4) = 


2S^S[JZ] 


(17) 


^^(2/1,222)^77(2/4,2/3) ' 

If S[TZ] is usual functional of the Hartree-Fock-Bogoliubov theory built up on the basis of 
the Hamiltonian which includes only two-particle density-independent interaction with the 
antisymmetrized amplitude we have: 

• 7 ^( 2 / 1 , 2 / 2 , — 2^Xx,X2^Xz2X4, 

x(^Ei,X 3 [ 4 i,+1-7^"^ (^ 11 ^ 2 ; x^.Xi) -f JF+(a; 2 ,a;i; 0 : 4 , 0 : 3 )] 

“4i,-X3 4xi,+i-7^^(3^1i^2; XA,x:i) + 5^^,-ir^{x2,xi] 0 : 3 , 3 ^ 4 )] ^ 

+ 4i,-X2 43,-X4 4i,X3 4x2,+i-7^“(^1i^ 2; )r3, 2 ) 4 ) + 5x2,-1-7^“(^3, 3 ) 4 ; , (18) 
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where 


T"^{Xi,X 2] = W^‘^\x2,X2,] Xi.Xi) , T [Xx,X2\ X^^.Xa) = [xx,X2\ X^.Xi^ . (19) 

In what follows we assume that Eq. m is fulhlled for the interaction JF, however Eqs. dH 
are not supposed to be fulhlled. In other words, the amplitudes JF+ and will be consid¬ 
ered to be independent, and it will be supposed that no other independent components are 
contained in the amplitude T. These assumptions correspond to the approximation adopted 
in the Theory of Finite Fermi Systems with pairing correlations (TFFSPC, Ref. jldji. 

3 TRANSFORMATION OF THE EQUATION 

FOR LINEAR RESPONSE MATRIX TO CHANNEL FORM 

The LRM dehned by Eq. ()12|1 contains information about excitations of the initial system 
in three different channels corresponding to the transitions to the states of the hnal systems 
with different numbers of particles. Suppose that the number of particles in the ground state 
of the initial system is conserved and is equal to Aq. Let A be the number of particles in the 
hnal system. Then, in accordance with the standard terminology, we have: (i) ph channel 
if A = Aq ; (ii) pp channel if A = Aq -|- 2 ; (iii) hh channel if A = Aq — 2. Notice that it is 
not necessary to introduce hp channel explicitly because of symmetry of the LRM and other 
quantities. 

Let us introduce the channel index c G {ph, pp, hh} and dehne the projection opera¬ 
tors 


[I] 

to 

CO 

“ <^X3,+i'^X4,+i ^3) < 5 (a; 2 , 0 : 4 ) , 

( 20 ) 

E^PP\xi,X2-, P3,P4) 

= ^X3,+i ^X4,-i (-1)^^'"" X3) S{x2, X 4 ) , 

( 21 ) 

E^^’^\xi,X2; P3,P4) 

= ^X3,-l^X4,+l X3)6{x 2, Xi). 

( 22 ) 


The sense of these operators is obvious: acting on any quasiparticle operator they cut out 
its components with hxed x values. Thus applying each of these projectors to LRM and 
correlated propagator we obtain the following components: 


a;) = 

j dys dyi E^^\xi, X 2 ] 2 / 3 , Vi) Hvs, Va, , 

(23) 

h}^^^\xi,X2] UJ) = 

j dy3dyiE^‘'\xi,X2; P3, P 4 ) A°(p3, P 4 ; a;), 

(24) 
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= j dy,dy,E^^\x^,X2; y,,ye) 

X [A{y5, ye; X3+,X4+; uj) - A{y5, ye; X4-,X3-; uj)] , (25) 


Ad^''P^\xx,X2; X3,X4; uS) 

= j dy5dyed''\xi,X2; y5,y6) 



X A(2/5, ye; X 3 +, X 4 -; uj) ^ 

(26) 


= j dy5dyed''\xi,X2; y5,y6) 



X My 5 , 2 / 6 ; ^ 3 -, X 4 +; uj) (-i)5+'^3 , 

(27) 


where the second channel indices in the Eqs. are hxed by the X"iiidices of the 

propagators on the right-hand sides of these equations. 

Let us also denote: 


T‘'''''''\xi,X2; X3,X4) 

= 5c,c' ^'^’'\xi,X2; X3,X4) + X3,X4), 

(28) 

:F^p’^\xi,X2; X3,X4) 

= J^"^{Xi,X2; X3,X4) , 

(29) 

^^PP\xi,X2; X3,X4) 

= (-1)5+'^2+H'^4 jf-(x 3,X4; a;i,X 2 ), 

(30) 

:F^^^\xi,X2; X3,X4) 

= (-l)5+‘^l + i+'^3 ^ 3 , 0 : 4 ) . 

(31) 


The additional restoring amplitude jg introduced in Eq. (EHl) for the purpose of 

“forced consistency” and will be specified in the following. 

Making use of the definitions (HHD, O-ijSID one can rewrite the Eq. (ED in the 
channel form 


A(")(a:i,X2; co) = k^^^\xi,X2; dx3dx4dx3dxe A^^'^'\xi,X2; xe.x^; u;) 

dc" ^ 

X d‘''''^'^{x3,X4;x3,x^d'^'^{x3,xe;u))i (32) 


where summation is performed over all the channels. 


4 EQUATION FOR PARTIAL COMPONENTS OF THE LRM 
IN CASE OF NEUTRAL EXCITATIONS 

In the present paper we solve the LRM equation dSl for neutral excitations. To separate 
the angular dependence in a spherically symmetric system we use decompositions in terms of 
spherical tensor operators. In particular, it is assumed that the operator of the renormalized 
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external field in the Eqs. (d, d, d has the form 


V^{xi,X2) = EjV(a:i,a;2) = 5(ri - rs) (n) , (33) 

LS 

where 

TjLSM{'n)a^a^ = (Lm5'/i| JM) YLm{n){aSf,)a^a^ , (34) 

mfi 

asfj. are the Pauli matrices in the tensor representation: 


= V2(2S + l)(-l)i-'. 


1 

2 


1 

2 


5 


(72 —(Ji 




For electric multipole excitations we use the standard ansatz: 


( 35 ) 


VSLsAr}^SjLSsar‘‘eW, ( 36 ) 

where is an effective charge in the center-of-mass reference frame: — Z (— 

= (1 - A-y + (Z - 1) {-A-Y- 

Let us denote in accordance with Eqs. d, d, d, and d: 

A^/^{xi,X2; uj) = - J dx3dx4A^''’P^\xi,X2; X3,Xi; a;) EjV(a;4, xg). (37) 

Let Aj],j be the solution of the Eq. (j^ with = Aj^\ In this case the partial compo¬ 
nents of the LRM and of the related quantities are dehned as 

= j dndxidx2T]^s^{n)^^^^6{xi,X2;r,T)A%{xi,X2;uj), (38) 

= dndxidx2T]^g^{n)^^^^5{xi,X2; r,T)A]^^{xi,X2; uj), (39) 


^L^SrL's'T'iz^A']^) = dryj du du' dxi dx2 dx^ dx4 6{xi, X 2 ; T, t) 6{x3, x^] v', t') 

X X3,X4; Uj) Tj l'S ', (40) 

where 

AlsmM.,., = (-l)"^**+"rris_„(n),,,,, (41) 

6 {xi,X 2 ]r,T) = 6 r^^rdr^^rd{ri-r) 6 {r 2 -r). (42) 
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Further, we assume that the effective interaction is determined by the following de¬ 
composition: 


X^,Xi) = 5{ri-r 2 ) 5{r^-r^) 


X 

X 


LSL'S'JM 


q~J (C,C') 

'^LSr^.L'S'rg 


(ri,r3) 




(43) 


where 


'T'J (^X') 
'^LSt, L'S't 


r{ph) 
•^S, tt' 

Tripp) 

•^S, rr' 


k=l 

Co (<550 [/ + (2<5r,r' “ 1) /^]+'^Sl [d + (2 ^ t , t ' 

X'itl = SsoK.r'Xi. 


1) s']) , 


(44) 

(45) 

(46) 


In Eq. (14411 the first term represents Landau-Migdal zero-range interaction of the 
TFFSPC. In the standard parametrization, the quantities Cq, g, and g' in Eq. (03) are 
constants. The functions /(r) and /'(r) are determined by the parameters f /j„, /g^., //„, 
and by the nuclear density in the ground state poir) by means of the ansatz: 


/(r) = fa + ( fin - fa) So('')/Po(0) , 


(47) 


and analogously for /'(r). For the interaction in the pp and hh channels we have: 


=C„/ln(c,/a. 


(48) 


where Cp is a constant, ^ is an energy cutoff parameter. 

The second term in the Eq. (P|) is a correction corresponding to the amplitude 
jrrest{c,c') (OHll . The similar correction arises in the calculational scheme with 

“forced consistency” which was developed in Ref. |T3] to obtain the spurious isoscalar 1“ 
state at zero energy in the non-self-consistent approach for the case when only ph channel 
is taken into account. Here the straightforward generalization of this scheme is presented 
in which the pp and hh channels are also included. Making use of the same method as in 
Ref. [H] we obtain: 


7 (c) 

5t, 1 


V) = 5a.r,h s. 


c,ph OSO 


dU^jr) 

dr 


7 (c) 
St, 2 


r] = 


J2ri‘‘!rAr)QtUr), 


(49) 
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( 50 ) 


" = 


/•LXJ 


cSt 


where 


^OO 

drW'‘‘A\^‘;jl,^(r,r'-,Qj„)F‘^'^^(r'), cJo^O. 


(51) 


c'5' 


In Eqs. gHD, U'^{r) is an anxiliary potential, the well depth of which is chosen to satisfy the 
condition: 

roo 

/ drr^[Z (r) - N (r) ] = 0 , (52) 

Jo 

which ensures the spurious state excitation probability to be equal to zero. 

Making use of the above dehnitions we can reduce the Eq. (IS 2 D in case of the spheri¬ 
cally symmetric system to the following equation for the partial components: 


- E 


VS'r’d L”S”d'd' 


poo poo 

/ dr'r'^ / 

'0 Jo 


dr" r"^ A 




where 


AjlItAA = 

AlSri'^) = ^c^phVjLSTir) 


/ dr' r' 2 Ai^Sr2'S'r A VjAs’r d) , 

L'S'c' 


^O(c') / /N 


(53) 

(54) 

(55) 


In terms of the partial components the Eq. (HEI) for the polarizability takes the form: 

poo 

n,(i.,) = (2 J + 1) 5^ / dr V°tl{r) A'jisrin ■ (56) 

LStc “'0 

It is important to note that the amplitudes and which dehne the single¬ 
quasiparticle basis functions according to Eq. have to be determined from the solution 
of the gap equation with the same interaction which enters the Eq. (jh.Sjl for the LRM 
through the Eqs. dm) and (gni). Namely, we have: 


“a = 


1 + 


E. 


Vx = sgn(AA) 


'1 


1 - 


^A 

E, 


(57) 


where is the eigenvalue of the single-particle Hamiltonian corresponding to the eigenfunc¬ 
tion (Pxd)y is the chemical potential. 


Ex = 



+ Al. 


fJjj 


(58) 
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The values of the energy gap Aa are determined within the BCS approximation from the 
equation: 


A. = -E 

(A') 


2 jA' + 1 Aa' 

dvr 2Ey ’ 


(59) 


where 


•^(AA') 


= 5. 


dr ( 


-^(A')( 




(60) 


5 CORRELATED PROPAGATOR 

IN THE COORDINATE REPRESENTATION 


Let us obtain an explicit formula for the partial components of the correlated propagator 
in terms of the reduced matrix elements. First of all, notice that the following relations 
are fulfilled for the matrix elements A 12 , 34 ( 01 ;) entering Eq. m in case of the spherically 
symmetric system: 


^ 12 , 34 (^A') — ^( 2 J + 1 ) Aj^2,34](‘^) 


JM 


X (-l)J 2-"^2 

^fl2,34] (“A^) = (2J + 1) ^ 12 , 34 ( 0 ;) 


( '■ 

J2 

J \ 

1 (-1)14-^4 1 

^3 

J4 

] 

\ mi 

-m 2 

M j 

1 1 

1 m 3 

—m4 

M / 


, ( 61 ) 


mim2m3m4 


X (- 1 ) 12-”^2 


(n 

J2 

J \ 

1 1 

^ .3 

44 


\ mi 

-m 2 

M j 

f 1 

1 m 3 

—7714 

M j 


, (62) 


where 1 = [ 1 ] = {(l), 77 i}, ( 1 ) = (Ai) = {ri, ni,/i, ji}. Notice that with these 

abbreviated notations we have: . We will use 

the antisymmetric form of equations for the correlated propagator (see Ref. [Hj for details) 
which allows to reduce the dimensions of matrices entering these equations. In this case the 
following relations are fulfilled: 

^fl2,34](^) = -^li^2(-l)^^''“'^^''”'^f2l,34](^) 

= -^3^4(-l)^^'^“'"^'^'^fl2,43](^)> (63) 


where [!] = {(!), -r]i}. 

Making use of Eqs. ®, dHD, (0, and m it is easy to show that the functions 'ipi{y) 
entering Eq. (d obey the equalities: 

^i{x]+) =w^ip^^{x), = r]^Wjip^^{x), (64) 
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where 


Wi = w^[i] = ,+i W(i) + ,-i ^^(1) • (65) 

Let us define the reduced matrix elements of the operator Tjlsm by the relation: 


{jihmi\TjLSM\j2hm2) — {—ly^ '"i | | {jili\\TjLs\\j2h) ■ 

nil —fn2 —M 

In the explicit form we have: 

{3M\TjLs\\32l2) = i [l + (-l)^+h+h] (_l)S+. 2 -^ 


X 


X 


(2J + 1) (2L + 1) (2ji + 1) (2j2 I 

to ' i -i 0 


J L S \ S{S + 1) I J L S 

0 0 0 /^V'^('^+l)\l 0 -1 


( 66 ) 


X [ (^1 “ ji) (2ji + 1) + (— (I 2 — j2) (2j2 + 1) ] f • (67) 


Making use of Eqs. (0, (0), (HH), (I2ni)-(|221), (I2SI)-(EZI), dED), dESl), (EH), dEEl), we 
obtain from the Eq. ( 001 ) the following ansatz: 


where 


^LSt,L' 5 't'(L ^ ) ^) ( ^ ^c',ph ) ^r,T' ^T^,T ^T^,T ^ 7 - 3 , r' ^(21) ^(43) 

[1234] 




iim(r) 

x(ofl2] - (-l)'’a[2l]) 4i'i2,34|(x') 

x(oj34] - %')4 (-Ij'^'afai). 


( 68 ) 


(69) 


T, 


( J) LS,L'S' 


(12,34) 


2J+ 1 


( 32h 11 Tjls 11 ji^i) (jih 11 Tjl>s' 11 jsh ), 


(70) 


“fl2] = ^2 ’ "fl2] = ^2 ^1 ^2 ’ "[12] = W 2 . (71) 

Theoretically, summation in the Eq. EHl) is supposed to be over complete ordered set 
of the states forming the doubled configuration space. This summation is facilitated due to 
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the symmetry defined by the Eq. (P|). The order-bounding factors are defined as follows: 
0 ( 21 ) = 1 if the ordinal number of the state ( 1 ) is lesser than the number of ( 2 ) [( 1 ) < ( 2 )], 
= i if (1) = (2), = 0 if (1) > (2). 

However, in our calculational scheme, summation in the Eq. dEHl) is restricted by 
the discrete and quasidiscrete states entering a valence zone near the Fermi level. It is 
supposed that the remaining part of the sum can be approximated fairly well by the RPA-like 
propagator which contains transitions from the quasiparticle levels to the continuum. 

Thus, we use the scheme which is analogous to the ones described in Refs. iHnni. According 
to this scheme the total correlated propagator is represented as a sum of two terms: 


aJ (c.c') 

^LSt,L'S't' 


r'; uj) = A 


J {c,c') cont 
LSt.US’t' 


r'] a;) -I- A 


J ic.c') disc / / 

LSt,L'S’t'\' y ' y 


U! 


(72) 


where 


w) = -5c,c'5c,ph5ry\ ^( 2 ) W ^(2)(0 

ii^i (2) 

-i(J) LS^US' V ^ nor 


disc 


X r," [G";y(r,r'; fi, - + w) + (-1)^+® G”')y(r,r'; fi, - B,.,, - x,) ] 


_ _ \S+S' /^nor 

qi2,i2) 

disc 

^ ^ , T , T ^(21) “^(l) '^(2) (^) “^(l) ) '^(2) ^ 

(12) 


xT, 


(J) LS,L'S' 


( 12 , 12 ) 


d 2 ) 


fil) 


.UJ + fir — -E(2) — £(i) (J — fJ,T + -E(i) -1- £(2) 


+ (- 1 ) 


S+S' 


hi) 


d 2 ) 


U) flj- — -^( 1 ) — ^[2) ^ ^(2) ^ 


( 2 ) ^^( 1 ) 


(73) 


disc 


ASA!'s‘"(’-y; y = ic.c-y,yr.r''£Sr,.rSr,.rffl’„ |l - O,),,] 


( 12 ) 


X R(i)(r) R( 2 )(r) R(i)(r') R( 2 )(r') [up) (1 - nfa)) + (1 “ ^^(p)] 


X 


(- 1 ) 


S+S' 


.u + E, 


( 1 ) 


E, 


( 2 ) 


-^( 1 ) -^( 2 ) 


+ (1 + ^c',p/7) ^r,T 


wind 

E 




J {c,c')LS,L'S', 


(21) ^^(43) ^[12,34] 


o; 


(74) 


[1234] 

In Eq. (ESI), G^i°^{r,r'] e) is the partial component of the normal single-particle GF without 
pairing. It is calculated via the regular and irregular solutions of the Schrodinger equation 
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that allows to take into account continuum completely on the RPA level (see, e.g., Ref. na). 
The relation between the Eqs. and the initial Eqs. (EHl), (El is determined by the 

formal decomposition: 


( 1 ) 




(75) 


e — £ 


( 1 ) 


In Eq. flTHl . are the bounding factors dehned as: = 1 if (1) < (A;,), = 0 

if (1) ^ (Afe), where Ab is the bottom level of the valence zone which coincides in our 
calculations with “pairing window”, i.e. with zone where the energy gap is not equal to 
zero (see Sec. |H1 for details). The superscript “disc” in the ^ symbols means summation 
over all discrete and quasidiscrete states of the single-particle spectrum. The superscript 
“wind” means summation over the discrete and quasidiscrete states only inside the “pairing 
window”, i.e. the same summation as in the Eq. El- 


6 THE AMPLITUDE OF THE QUASIPARTICLE-PHONON 
INTERACTION AND ITS MATRIX ELEMENTS 

Matrix elements of the amplitude of quasiparticle-phonon interaction belong to the basic 
quantities entering formulas for the correlated propagator of the QTBA (see Ref. [0] and 
Sec. 13 of the present paper). In accordance with the general dehnition of this amplitude in 
the microscopic approach we have: 

g'^{yi,y2) = j dy3dy4J^iyi,y2; y3,y4) p""^{y3,yi), (76) 

where the interaction amplitude JF is dehned by the Eq. El, is the transition amplitude: 


p"*°( 2 /i, y2) = {m\ b\yi) 6 ( 7 / 2 ) | 0 ). (77) 

The operators b{y) = b{x, x) are dehned by the relations: 

b{x,+)=a{x), b{x,—)=a\x), (78) 

where al(x) and a{x) are creation and annihilation operators of particles. 

As a hrst step we consider the coupling of the quasiparticles with core excitations only 
in the ph channel. So, substituting Eq. El into Eq. (El, we keep only the contributions 
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of the amplitude Taking into account the symmetry: 2 / 2 ) = —l/i), we 

obtain from Eqs. m and (iznD: 

^™(2/i, 2/2) = ^Xi,X2 , (79) 


where 


g"^{xi,X2) 


j dx3dx4T'^{xi,X2; 0 : 3 , X 4 ) p™°(a; 3 +, 0 : 4 +). 


In what follows we assume that the function g"^{xi,X 2 ) has the form: 


(80) 


j”(ii,X 2 ) =9“"( 2 ^ 1 .^ 2 ) = S{ri - r 2 )y,r, 5^J?y(ri) U,LSM(ni),ja, , (81) 

LS 

where we separate the projection of the total angular momentum M from the remaining 
phonon quantum numbers, which are denoted by the index q (including the total angular 
momentum Jq). Notice that Eq. CD is in accordance with Eq. ini) and with the decompo¬ 
sition of the effective interaction within the TFFSPC 
Let us denote 


= J dyidy2^i{yi)^2iy2)g'^^{yi,y2), 

9l^x^ = j dxidx2^),^{xi) ^1^{X2) g‘>^{xi,X2). 

From Eqs. (0, ii, iH), and il it follows that 


9?" = (-!)*-”« 


J2 Jl Jq 
m2 —mi —M 



where 


(82) 

(83) 


(84) 


dm = E) 9('i“{92<2|IU,is||iili), (85) 

L 

poo 

9112) = R{i)ir) %)ir) gl^^{r) ■ (86) 

V 0 

Making use of Eqs. dm, ii, ii-ii, and the symmetry of the operator T jlsm we obtain 
the following two main formulas for the matrix elements of the amplitude of quasiparticle- 
phonon interaction: 


9I2 = 


32 jl Jq \ q 
m 2 —mi —M j 


(87) 
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where 


( 88 ) 


®fl2| (-1)® “"2 ) 9(f2) ■ 

s 

The Eq. ()88|1 can be easily generalized to take into account the coupling of the quasi¬ 
particles with core excitations in the pp and the hh channels. Proceeding in the same way, 
we obtain: 


^?i2] = ("[12] “ ^c,ph V 1 V 2 (-l)^«[2i]) {J2h 11 Tj^ls 11 Jih ) > (§9) 

LSc 


where the amplitudes defined by the Eqs. (f7T|l . 


/ dr (r) (r) (r) . (90) 

0 

Notice that in the Eqs. (IHH) and (jHHjl we have: To determine the 

function gT^^^^\r) entering Eq. (inni) the following method is possible. In terms of the partial 
components the Eq. dZSD reads: 


gqLSic) (r) (r) , (91) 

t' 

where functions are determined by the Eqs. (gSl) and is the partial 

component of the transition amplitude, and the contribution of the restoring amplitude in 
Eq. (11^ is omitted. From the spectral decomposition of the partial component of the LRM 
one can End that up to irrelevant constant phase factor the following equality is fulhlled: 




B{q)] J 2m 

Ca 


(92) 


where ^ solution of Eq. B{q)'\ = B{g.s. —>• q) is reduced excitation probability 

which is determined by the formula 

c. 

The polarizability nj(a;) in Eq. (innD is determined by Eq. dSEI)- In the Eqs. (P|) and the 
integration is performed along the contour Cg in the complex energy plane which encloses 
the poles of the functions Aj^^^^(r; a;) and nj^(a;) corresponding to the phonon energy 
ujq. This method automatically ensures correct normalization of the functions and 

^qLS{c)^^y However, to determine these functions in case of the excitations in the pp and the 
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hh channels one has to solve Eq. (HI for the different external helds Vj^g^{r) which take 
non-zero valnes for c = pp and c = hh in contrast to the definition dSSI). 

Notice that in order to calcnlate the fnnction gT^^^^\r) by means of the formnlas 
dnu-dni it is not necessary to know the exact valne of the phonon energy. It is snfficient 
that the pole uj = ujg (and only this pole) wonld be inside the contour Cq. The exact value 
of ujq can be determined by the same method with making use of the formula 

1 f du ^ 

^ Wh J ’ 

c. 

Another method to determine matrix elements consists in solving QRPA equation 
for the transition amplitudes making use of the basis restricted by discrete and quasidiscrete 
single-particle states entering a valence zone near the Fermi level. In terms of the reduced 
matrix elements this equation reads: 

wind 

i^LOq — p [E(i) -|- -E(2)])P(i2)j? ~ ^*^(12)??, ( 34 ) 77 ' ^(43) d(34)?7' ’ (9^) 

77 ' (34) 

where 


rJ 

•^(12)77, (34)77' 


[12,34] 


= y s s 


1 


^%-V' -^[12,34] 7 


( 96 ) 


("[12] “ ^1^2 (-1)^«[2I]) («[34] - V 3 V 4 (-1)^«[43]) 

LSc 

vfl-UA \ rpiJ) LS,LS q:.ic)S 
^ ^ ^c,ph ) (12,34.) '^(12,34)7 


( 97 ) 




■{c)S 


(12,34) 


= h 


'ri,r 2 0r3,T4 


drr‘^R^^^{r) . ( 98 ) 

The quantities entering these formulas are defined by the Eqs. (03)7 (Iini)7 dZi, dUD- Notice 
that QRPA equation (I95|l includes contributions of all the channels: ph, pp, and hh. 
Normalization condition for the transition amplitudes P(i 2 ),, has the form 


wind 

EE ^(21) ^ |Pp2)77l 2Jq + l. 


( 99 ) 


V (12) 


The reduced matrix elements and the reduced probability B{q) ] of excitation induced 
by the external field dMD are determined by the found values of p\i 2 )r^ through the formulas: 


wind 


9[12] ^2 •^[12,34] ^(43) , -V 4 ^'( 34)773 7 

[34] 


( 100 ) 
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B{q)] 


1 

2Jq + 1 


wind 2 

X] X] ^( 21 ) i^Jg) (21)71 ^( 12 )»? ’ 

?? ( 12 ) 


( 101 ) 


where 

i^j)(21)n = 5Z^^[“(1)^(2) + (-1)^^(1)“(2)] (^2/2||Tjls||jJi)(^JLs)(21)> (102) 

LS 

poo 

i^J73)a2) = £irr^fi,„(r)fip,(r)l>,V^(r). (103) 

d 0 

7 CORRELATED PROPAGATOR OF THE QTBA IN TERMS 
OF THE REDUCED MATRIX ELEMENTS 


In this section we draw the formulas for the reduced matrix elements of the correlated 
propagator Aj ^2 34 ]('^) which are obtained within the QTBA in case of the QPC model. In 
what follows, the summation over single-quasiparticle indices means summation over the 
discrete states inside the “pairing window”. The general formula for the propagator of the 
QTBA satisfying Eq. reads: 


^[12,34] (^) “ [*^[15] '^[26] + Q[12,56] (^) ] ^[5l78/(^) [*^[73] '^[84] + ^[7l34/ (‘^) ] 


[5678] 




(104) 


where 


= «(43) [QSV)+ (-!)"+'»-'• (!*>)]. (106) 

«(21, [Q|iUP) + QptalP) ] • (1“) 




—)o.O — A A A -^(—) 


[12,34] 


( uj ) - 6 6 A 


Vi -V2 "iJs -’?4 (12)i?i, (34 )i?3 


(^) 


The 


quantity \ 34 )r]' ('^) ^ solution of the equation 


^(12)i?, (34 )j 7' (^) ^(12)r), (34)r)' (^) +EE ^(65) ^(I2)v, (56)1?" 


(o;) A 


j(—) 
(56)»?", (34)?) 


(107) 


(108) 


T)" (56) 


where 


A 


j 

(12)71, (34)7?' 


i(^) = 


V^TI' [<^(13) <^(24) + (-1)^+'^ S(M) <^(23)] 

2 {u - T] [E^i) + E^2)]) 


(109) 
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^{12)»?, (34)r?' (^) ~ 


UJ- 7] [E(i) + E(2)] 


( 110 ) 


*^(12)7?, (34)7?'(^) X] '^7?j,r?'^7?2 -7?'^»?3,r?''^7?4 -7?'*^[12,34](^) • (m) 

It is easy to see that solution of the Eq. (HDHl) possesses the following symmetry: 

= (112) 

It enables one to determine all the elements of the matrix )( 34 )r,'(^) by solving Eq. (HUHl) 
for the non-zero block of the matrix 6 '( 2 i) ^( 12 ) 7 ? ( 34 )r?'(‘^) ^( 43 )- However, as follows from Eqs. 
dnnD, (dD, (HED-diiiii), for the construction of propagator in the coordinate representation 
it is sufficient to determine the elements only of this block. 

In order to dehne the remaining quantities which enter Eqs. (Uni-dinsi), (cni) in 
case of the quasiparticle-phonon coupling model let us introduce notations: 


Hi 2 g — -^12 + hi 7 -^12 ~ hi [ -^( 1 ) + -^( 2 ) ] 7 


(113) 


n'J =A n^* -L A i — l y'l+t2+t3+l4 77'^* 

-^[12,34] 7? h[i3] y[24] + '^7?,-! I -‘-J h[3i] y[42] 7 


(114) 


X, 


Jq _ 

[12,34] 77 


— ^- '\^J+Jq+j2 H 


3i 32 J 


D 


7 i 

Y,t .7 = dL _ , (115) 


f ""[12,34]7? ’ "[12,3] 2?' -H [12,33]7?3 

34 33 Jq \ ^3i + ^ 


where the reduced matrix elements are dehned by the formulas of Sec. IHl Making 
of the Eqs. pi3j) - ()115jl we obtain: 


use 


$ 


[12,34]' 


(-)=«>f.S(-)++I’l'sv). 


jJ(GSC)^ 

qi2,34] 


(comp), 
[12,34] 


(116) 


jJ(res) 

■ [12,34]' 




hi ^»?3 -'?4 


<? 




'?l -^3 VcU-fl 


J q W^ 

[12,34] 7?3 [12,34] 7?2 


32 q 


00 — Vt 


14 q 


^ b^[13,5] , ^ '^’?6-^2 b^[42,6] 

^[24] +^[13]Z^ 


[5] 


52 (J 


[ 6 ] 


16 (J 


( 117 ) 
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® [12,34] 




yJ <? Y'^ 

. [12,34] ,7^ [12,34] r,2 

I ^42, 


‘13g 


+ ^[24] 

[5] 

+ '^[13] Y 


^[13,5] ~ -^12 ~ ^35 q) 

^15q ^35g 


^Vfi,Vi ^[L,6] -^12 ^64(j) 

^62 q ^64 q 


( 118 ) 


[56] gg' 


^%<V2 \’Vi ^[13,5] ^[L,6] 
^15 q ^35 q ^62 q' ^64 g' 


X (a; + ii^56 — ^12 q — ^34 q ') , 


(119) 


pj(++), . 

-^[12,34] 


^1 '^^ l .»?2 


r A 1 

I % ->?3 [12,34] 

^13 g ^24 g 


U — il. 


32 q 


U — il 


14 g 


'^[ 24 ] Y 


%-Vi ■'[13,5] 


[5] 


^51 g ^53 q{^ ~ ^52 g) 


'^[ 13 ] Y 


% .-’?2 [42,6] 


[ 6 ] 


f^26 g ^46 g (itl — ^16 g) J 


, ( 120 ) 




= 5. 


A ^ A X'^ ^ 

'Jl --^3 [12,34] »73 [12,34] »73 


’?l,»?2 '?3’ ’?4 (a; — f232q) f224g (l^ ~ ^14g) 


5 


5 


[24] Y 

[5] ^ 

[13] Y 

[ 6 ] 


A A 

"’Is-y ■'[13,5] ^ "’?5,-»?i ■'^[13,5] 


-£'13 ^53 g 




51 g 




+ 


13 


A y^ A y? 

■'[42,6] _J_ "»?6,-»?2 ■'[42,6] 


-E 42 ^46 g 




26 g 




+ 


U — ^52 q 
1 


42 


U — fl 


16 g 



, ( 121 ) 


y (—^)/ 

' [12,34] 


QflO Q/,1 {^ ) 




Ki'-V2 


n ^ 


+ ^[24] Y 
[5] ^ 

+ ^[13]Y 


A X'^ "l A X'^ 1 

'^»?4--’?2 [12,34] r;^ ^ '^»?3 ,-^1 ^ [12,34] rjj 

^42g ( 1 ^ — ^32g) ^13 g (l^ ~ ^14 g) 

1 1 


^[13,5] ^ ^[13,5] . _ ^ 

-£'31^51g ^53 g \ -£'31 — 1152 g 


‘53 g 


[ 6 ] ^ 


A y^ A 

^^6 '^4 -'[42,6] ^Vq-Vj [42,6] . ^ 

-£^24 II 26 g II 46 g \ -£^24 ^ — ^36 q 


‘46 g 



( 122 ) 
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The Eqs. ()104|1 - ()122|) completely define the reduced matrix elements 34 ] (<^) which 
enter Eqs. (innD and (m for the “discrete” part of the total correlated propagator in the 
coordinate representation. Notice that in the Eq. (IHn]) the matrix elements 34 ] (<^) cire the 
same for all the channels which differ from each other only by the amplitudes (A^i 2 ]- Notice 
also that the form of the above equations for the reduced matrix elements y4j^2 34] ('^) 
same both for magic and for open-shell nuclei. In the former case, however, the index 7]^^ 
in the set [ 1 ] = is not independent, but is determined by the occupation number 

= 0 or 1 of the state ( 1 ) = (Aj^) as: rj^ = 1 — 

8 CALCULATIONS: DETAILS AND DISCUSSION 

As an application of our approach isovector electric dipole excitations in the GDR region 
have been calculated in semi-magic isotopes Ii6,i20,i24gj2. We started from a description 
of the independent single-particle motion in the standard phenomenological Woods-Saxon 
potential. Then the gap equation dSHD was solved for neutron subsystem under the usual 
condition that the number of particles in valence zone is conserved on average. In our 
calculations for the tin isotopes this zone (“pairing window”) consists of all discrete (^a < 0 ) 
and quasidiscrete states above the chemical potential /i and of the discrete states below 
/i starting from I/7/2 subshell for neutrons and 1 ^ 5/2 subshell for protons. We emphasize 
that the same valence zone restricts summations both in the gap equation dSHl) and in the 
equations of Sec. [3 for the QTBA propagator. 

As the quasidiscrete states we chose the discrete states with £a > 0 calculated with 
a box boundary condition and having the last extremum of the radial wave function in the 
range of non-vanishing values of the discrete spectrum (^a < 0) wave functions. According 
to this criterion the quasidiscrete states are selected quite well if size of the box is not too 
large. Notice that the partial components of the normal single-particle GF r'; e) 

entering Eq. dZSD were calculated with outgoing wave boundary condition except for the 
small number of components for which j7-values coincide with j7-values of the quasidiscrete 
states. For such partial components the box boundary condition was used in order to avoid 
inconsistency in the calculations. 

The parameter Cp in Eq. was chosen so as to obtain the averaged solution of the 
gap equation Aa to be equal to the averaged empirical value (see Ref. [in]): 

A = 12 MeV X A -^/2 _ (^^23) 
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For we have obtained A 1.1 MeV with Cp = 0.719 MeV. This parameter was used 
for the remaining two nuclei. The quantity ^ in Eq. is determined as follows: ^ = a/A ^2 
where A = h ~ '^2 = max(e;^) — /i. 

We assumed that the observable energies of single-particle excitations in the neigh¬ 
bouring odd nuclei have to be reproduced in the framework of the mean held plus BCS 
model. These observable energies (experimental energy differences) are dehned by the fol¬ 
lowing equation: = ± [T;^(A ± 1) — T(^g(A)] , where is the ground state energy of 

the even-even nucleus (core) consisting of the A nucleons, is the energy of the ground 
or the excited state of the neighbouring odd nucleus consisting of the Ail nucleons. In 
order to get an agreement with the experimental energies the well depth of central part of 
the Woods-Saxon potential was slightly varied so as to obtain for protons and 

i Ex for neutrons where is connected with by the Eq. (ISHD. Thus obtained 
energies and wave functions form the above mentioned basis set </3 a(3^)}- 

In the present calculations we included ground state correlations (GSC) only in the 
QRPA part of the correlated propagator. Another type of GSC caused by quasiparticle- 
phonon coupling (GSC/QPC) and originated from backward-going terms in time-ordered 
diagrams was not incorporated. This means that: (i) only the term (a;) in the 

Eq. flll6|) is accounted for; (ii) associated components of the propagator containing the 
functions (<^) [see Eqs. p04|) -p06j). p2()j) -p22|)] are 

excluded from the calculations. 

Since our single-particle and single-quasiparticle energies are htted to the experiment 
as it was mentioned above, these energies and corresponding wave functions already contain 
effectively the admixture of phonons. This phonon contribution should be removed from 
the mean held, energy gap, and the effective interaction to avoid double counting if the 
quasiparticle-phonon coupling is included explicitly. To solve this problem we use the method 
which corresponds to the self-consistent scheme of the QTBA (see Ref. jH]) if the GSC/QPC 
are not included. In this case the method consists of the replacement of the amplitude 
<hj ^2 34 ](‘^) the Eq. (Illl|l by the difference amplitude |)j^ 2 '^ 34 j(a;) where 

4i'SH = 4S(-)-«i'S(0)- (124) 

Making use of this subtraction procedure, ground-state contributions of the quasiparticle- 
phonon coupling dehned by the quantity d>j^ 2 '^ 34 j( 0 ) are removed from both the mass operator 
and the ehective interaction. 
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In the calculations of the GDR we neglected the contribution of pp and hh channels 
that is justihed for this type of excitations (see Ref. [IH])- The parameters of the Landau- 
Migdal zero-range interaction in the ph channel entering Eqs. (Hi and were taken in 
accordance with the standard set which is usually used in similar calculations (see, e.g., Refs. 

mm)- 

/,„ =-0.002, /L = 2.62, 4 = 0.76, 

g = 0.05, g' = 0.96, Cq = 300 MeV fm^. (125) 

Notice, however, that the parameter 4 has been changed as compared to one used in 
Refs. ^ E] in order to reproduce the experimental mean energy of the GDR in ^^°Sn 
within QTBA calculation. The parameter 4 was slightly varied from —2.187 to —1.957 in 
the computation of the energies and amplitudes of low-lying collective 2~^ and 3“ phonons 
within QRPA in conhguration space. This allowed us to £t the energies of these phonons to 
experimental values. Then the value of fex was averaged and used in the calculations of the 
remaining phonons and further in the QTBA calculations. 

The phonon characteristics were calculated within the QRPA by making use of Eqs. 
(PB -lfTUH)). The conhguration space was restricted by the “pairing window” dehned above. 

(c)S 

In this calculation the matrix elements of the interaction 412 34) with S = 1 were omitted 
since they do not give signihcant contribution in case of low-lying collective modes with 
natural parity. In addition, we have neglected contribution of pp and hh channels. These 
simplihcations decrease dimension of the QRPA matrix in the configuration space by a factor 
of two. For all three chosen tin isotopes the collective modes with spin and parity 2~^, 3“, 4~^, 
5“, b’*', and with energies below the neutron separation energy were included into the phonon 
space. A mode is assumed to enter the phonon space if its reduced transition probability 
is more than 10 % of the maximal one for hxed spin and parity. This value was taken as 
an approximate criterion for phonons selection. Characteristics of the low-lying vibrations 
taken into account in the calculations are collected in the Table ^ 

The dipole photoabsorption cross section aEi{E) is the basic observable computed in 
the present work. This value is expressed via the strength function Sei{E) according to the 
well known formula: 

aEiiE) = -^ESEi{E). (126) 

'dric 

The strength function in turn is simply connected with the polarizability n£;i(a;) [see Eq. 
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Table 1: Energies and reduced transition probabilities of the low-lying phonons used in 
QTBA calculations. 


J" 

ii^Sn 

u 

(MeV) 

B(EL)t 

(e^fm^-^) 

i20Sn 

u 

(MeV) 

B(EL)t 

(e^fm^-^) 

®24Sn 

u 

(MeV) 

B(EL)t 

(e^fm^^) 

2 + 

1.29 

2.97x10^ 

1.17 

3.53x10® 

1.13 

2.55x10® 

3- 

2.27 

9.06x10^ 

2.40 

8.42x10^ 

2.60 

1.32x10® 


5.21 

1.35x10^ 

5.07 

1.89x10^ 



4+ 

2.60 

1.59x10® 

2.62 

1 . 22 x 10 ® 

2.11 

1.19x10® 


6.56 

5.15x10® 

7.58 

8 . 12 x 10 ® 

3.93 

5.87x10® 


4.88 

3.08x10® 

5.35 

3.35x10® 

7.60 

5.33x10® 


5.40 

3.02x10® 

4.89 

3.26x10® 

5.06 

4.17x10® 


5.56 

2.58x10® 

3.44 

2.25x10® 

7.36 

2.31x10® 






5.38 

1.64x10® 






5.59 

1.35x10® 






4.48 

1.30x10® 

5“ 

5.91 

2.95x10^ 

2.67 

1.98x10^ 

6.85 

2.85x10^ 


2.75 

1.46x10^ 

7.04 

1.80x10^ 

2.40 

1.47x10^ 


6.13 

6.53x10® 

5.52 

1.41x10^ 

5.49 

1.41x10^ 


3.49 

5.52x10® 

6.53 

9.98x10® 

4.03 

1.14x10^ 


8.39 

4.15x10® 

3.56 

8.57x10® 

8.48 

9.66x 10® 


8.98 

3.57x10® 

8.59 

8.43x10® 

6.09 

5.05x10® 


3.82 

3.13x10® 



4.38 

4.95x10® 

6 + 

5.38 

7.63x10® 

7.44 

9.74x10® 

4.78 

1.04x10® 


5.70 

7.51x10® 

5.34 

7.56x10® 

7.24 

8.79x10® 


4.84 

5.93x10® 

4.27 

5.63x10® 

5.06 

5.59x10® 


6.41 

5.48x10® 

4.88 

4.33x10® 

2.40 

5.47x10® 


3.35 

4.29x10® 

2.94 

3.93x10® 

5.31 

5.46x10® 


3.99 

2.27x10® 

7.83 

1.80x10® 

7.67 

2.58x10® 


6.74 

1.32x10® 

4.98 

1.67x10® 




3.70 

1 . 00 x 10 ® 

4.44 

1.58x10® 
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Sei{E) = -hmUEi{E + zA). (127) 

TT 

As follows from the Eq. dsni), to determine the value of function n£;i(a;) at a given complex 
energy variable oj and, consequently, to compute the values of Sei{E) and aEi{E) one has to 
solve Eq. for the LRM. This equation was solved within the framework of two models: 
QTBA and QRPA. The spurious isoscalar 1“ state has been eliminated using the “forced 
consistency” method presented in Sec. El The single-particle continuum was included as 
described in Sec. El Following calculational scheme which is usually used in the response 
function formalism, in the present work we introduced a smearing parameter A which is 
actually an imaginary part of the energy variable oj. In the calculations of GDR the value 
A = 250 keV was used. 



Figure 1: El photoabsorption cross section for ^^®Sn calculated within QRPA (dashed line) 
and QTBA (solid line). The smearing parameter A is equal to 250 keV. 
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Figure 2: Same as Fig. [H but for 


Calculated pliotoabsorption cross sections for three above indicated tin isotopes are 
shown in Figs. ^0 The solid and the dashed curves represent the QTBA and the QRPA 
correspondingly. In order to compare the calculated cross sections with experimental data 
the parameters of Lorentz function <Jl{E) were found where 

F2^2 

(E2 - ElY + F2E2 ■ 




The parameters Eq, F, and ag in the Eq. (I128|l were obtained by making use of the following 
condition: the energy-weighted moments mg, m_i, and m _2 of the functions (insD and 
p28j) should coincide. This method is analogous to the one developed in Ref. UTI but in 
contrast to Ref. UH we calculated the moments in the hnite energy interval 0-30 MeV for 
which experimental data are available. Parameters of the Lorentz £t are compiled in the 
Table El In this table values of the depletion of the energy weighted sum rule (EWSR), i.e. 
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E [MeV] 


Figure 3: Same as Fig. [H but for 

integrated cross sections aEi(-E) dE, are also presented as percentages with respect 

to the corresponding Thomas-Reiche-Kuhn (TRK) values 59.74 NZ/A MeV mb. 

In the Table El experimental Lorentz function parameters of the GDR are shown. 
Experimental mean energies demonstrate the well known property to decrease against neu¬ 
tron excess. Our QTBA results drawn in the Table |21 reproduce quite well these mean 
energies. Since we take into account a hnite number of the low-lying phonons, our theo¬ 
retical curves have the shapes which are rather far from single Lorentzians. Nevertheless, 
remaining parameters of the Lorentz function calculated within QTBA, i.e. F and uo, are in 
a reasonable agreement with experimental values. 

Consider the results obtained within QTBA and QRPA to analyze effect of the QPC 
on the integral characteristics of resonances. As one can see from the Table El QTBA gives 
signihcant increase of the total width as compared to QRPA ( Fq-pg^ > 1.6 F^^p^) owing to 
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Table 2: Lorentz function parameters of the GDR in ii6,i2o,i24gj^ obtained within two 
microscopic approaches for 0-30 MeV energy interval. Values of the depletion of the EWSR 
are presented as percentages with respect to the corresponding TRK values. 



QRPA 

QTBA 

QRPA 

i20Sn 

QTBA 

QRPA 

QTBA 

Eo (MeV) 

14.74 

15.44 

14.65 

15.39 

14.35 

15.10 

P (MeV) 

2.4 

4.0 

2.7 

4.4 

2.6 

4.4 

(To (mb) 

452 

302 

423 

288 

452 

298 

EWSR (%) 

94 

99 

97 

102 

97 

102 

Table 3: Experimental Lorentz function parameters of the GDR. The results 

Refs. [IHI and [TU]. 

are taken from 



i20Sn 



ng 

PI 

na 

PI 

PI 

PI 

Eo (MeV) 

15.57 ±0.1 

15.67T0.04 

15.38T0.1 

15.40T0.04 

15.29±0.1 

15.18T0.04 

P (MeV) 

5.21±0.1 

4.19T0.06 

5.25±0.1 

4.88T0.06 

4.96±0.1 

4.81±0.06 

cTo (mb) 

270±5 

266±7 

284±5 

280±8 

275±5 

283±8 


contribution of the spreading width Th Clearly this result could be expected from physical 
point of view. The EWSR values obtained within the QTBA for the investigated energy 
interval 0-30 MeV are rather close to the TRK ones but again there is a noticeable difference 
between the QTBA and the QRPA values. However this difference has another source which 
is the subtraction procedure described above [see Eq. (I124|l ]. It can be rigorously proved that 
for the version of QTBA in which GSC/QPC are not taken into account and the subtraction 
procedure is not applied the equality EWSR^^pg^ = EWSR^^^^p^ is fulhlled exactly. Notice 
that the analogous equality is fulhlled between the values of EWSR dehned within the 
MCDD (Ref. HU]) and within the RPA. In our calculations just the subtraction procedure 
determined by the Eq. (ini gives rise to increment of EWSR in the QTBA. Switching off 
this procedure we have obtained that the values of EWSR calculated within the QTBA in 
the “inhnite” energy interval (0-200 MeV in our calculations) are equal to the corresponding 
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QRPA values with sufficiently high accuracy. This result can be considered as a test of our 
calculational scheme. The mean energy shift of about 0.7-0.8 MeV obtained in QTBA with 
respect to QRPA has the similar nature: it arises mainly due to the subtraction procedure. 
Without this subtraction the QTBA mean energies decrease as compared to the QRPA ones 
by about of 0.2 MeV. Thus, the subtraction procedure results in the signihcant change of the 
averaged characteristics of the excited states calculated in the QTBA as compared to the 
QRPA. On the other hand, it ensures elimination of the QPC contributions from the ground- 
state characteristics leading to the equality between QTBA and QRPA response functions 
at a; = 0. 

Finally, notice that although we have taken into account the most important effects 
of the QPC, the neglected contributions of the GSC/QPC and of the dynamical coupling 
to the pp and hh channels may also affect the results. In case of magic nuclei, the role of 
the GSC/QPC in the description of nuclear excitations was investigated in a series of papers 
(see, e.g.. Ref. [T] and references therein). The study of these effects in the nuclei with 
pairing is in progress. 

9 CONCLUSION 

The quasiparticle time blocking approximation (QTBA) is applied to describe El excitations 
in the even-even open-shell spherical nuclei. Within the QTBA pairing correlations, two- 
quasiparticle (2q), and 2q®phonon conhgurations are included. The model is based on 
the generalized Green function formalism in which the normal and the anomalous Green 
functions are treated in a unihed way in terms of the components of generalized Green 
functions in a doubled space. To determine response of the spherically symmetric system 
against an external held within the QTBA the integral equation for the partial components 
of the linear response matrix in the coordinate representation has been obtained including 
coupling between particle-hole, particle-particle, and hole-hole channels. Gonhgurations with 
a particle in the continuum are included into the QRPA part of the response function. This 
enables us to describe both spreading and escape widths of nuclear excited states. In our 
calculations we use phenomenological Woods-Saxon single-particle input and independently 
parametrized effective interaction of the Landau-Migdal form. So the additional procedure 
to eliminate spurious dipole mode has been formulated in terms of QTBA. 

The developed method has been applied to calculate isovector El strength distribution 
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in nuclei ii6,i2o,i24gj^ results for the photoabsorption cross sections in the indicated 

tin isotopes are presented. Since our main purpose was to test new approach, these first 
calculations have been performed assuming some additional simplihcations of the model: 
ground state correlations caused by quasiparticle-phonon coupling and dynamical coupling to 
pp and hh channels were ignored. Nevertheless noticeably fragmented giant dipole resonances 
have been obtained for all three investigated nuclei. Calculated integral characteristics of 
the resonances are in a reasonable agreement with experimental values. 
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